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INTRODUCTION 

In the theory of the strong interaction, Quantum Chro- 
niodynamics (QCD), it has so far been impossible to 
derive hadron masses directly using analytical methods. 
Lattice QCD is an approach whereby the functional in- 
tegral is evaluated numerically to obtain various physical 
quantities. One difficulty in this has been the inclusion of 
dynamical quark effects, but this is an area where much 
recent progress has been made. However, for computa- 
tional reasons, the masses of the valence quarks can be 
much more easily varied than those of the sea quarks. 
This situation is referred to as partially quenched lattice 
QCD. 

The quark masses which can be accessed by present 
simulations are, especially for the sea quarks, signif- 
icantly larger than the physical up and down quark 
masses. In order to obtain the masses of the lightest 
hadrons, one thus needs to extrapolate from the quark 
masses used in the lattice calculations to the physical 
ones. 

In this domain Chiral Perturbation Theory (xPT), 
which is an effective low-energy field theory approxi- 
mation to QCD, should be valid. It provides the cor- 
rect quark mass dependence of the pseudoscalar meson 
masses, including the nonanalytic dependences often re- 
ferred to as chiral logarithms. The three flavour xPT, as 
formulated by Gasser and Lcutwyler [jj], is valid in the 
QCD case with equal valence and sea quark masses, and 
cannot be used to extrapolate from partially quenched 
lattice QCD calculations. The extension of xPT to the 
quenched |2| and partially quenched Q cases has been 
carried out by Bernard and Golterman. The most general 
expression at one loop for the pseudoscalar meson masses 
was worked out by Sharpe and Shoresh [J|. In this work 
we present the result for the meson masses at next-to- 
next-to-leading order (NNLO) in PQxPT. It should be 
pointed out that the case of full QCD corresponds to the 
special case of partially quenched QCD where the sea and 



valence quark masses are put equal. The free parameters 
or low energy constants (LECs) of xPT, i.e. the QCD 
case, may thus be unambiguously determined from those 
of PQxPT. These points are discussed in detail in [J|. 

We work here in the version of PQxPT without the 
supersinglet $0, discussed in detail at one loop in p. 
Here we present the expression for the charged, or off- 
diagonal, meson mass for the case of equal valence and 
equal sea quark masses, but these different from each 
other. The charged meson masses for the more general 
quark mass cases, as well as more details of the present 
calculation will be presented elsewhere |a| . Planned work 
includes the decay constants and the neutral, or diagonal, 
meson masses. Note that the presence of NNLO terms is 
already seen in the works of ref . |3 . 

In the next sections we present the technical back- 
ground, the full expression for the pseudoscalar meson 
mass, and finally a short discussion and numerical results 
as a function of the input quark masses. 



TECHNICAL OVERVIEW 

Most of the technical aspects required for calculations 
at the two-loop level in PQxPT without $0 already ex- 
ist, but they have not been mentioned as such. The full 
divergence structure at one loop was worked out in ja, Isi] ■ 
There, it was noted that the divergence structure is re- 
ally identical to that of normal xPT with the number of 
flavours equal to the number of sea quark flavors, pro- 
vided that the traces and matrices are replaced by the 
supertraces and matrices appropriate for quenched xPT. 
This is justified by the replica method. In ref. [l(]j it has 
been argued that the quenched approximation, i.e. no sea 
quarks, can also be obtained using the replica method. 
This entails a calculation with np flavors, setting np = 
in the final answer. In ref. jlQ| this was proven explicitly 
at one-loop order. These arguments can be generalized 
to higher orders and to PQxPT. This allows the use of 



the known results for np flavours in normal xPT (va- 
lence and sea quark masses equal) , in order to obtain the 
needed Lagrangians at O(p^) and O(p^) from ref. [Il|, as 
well as the full divergence structure from ref. [iJI ■ In par- 
ticular, the U p fl avor Lagrangian at 0{p'^) contains the 
term with Lq [11|, see also [l|. In the context of PQxPT, 
this extra term was pointed out in ref. ^^. We will use 
the notation of refs. |llLll2l| for the Lagrangian at 0{p'^) 
and 0(/), L[,i = 0, . . . , 10 and iJ.[,i = 1,2 for the 
renormalized 0{p'^) LECs, and K^, i = 1, . . . , 115 for the 
renormalized LECs at 0{p^). The reason why these La- 
grangians can be obtained by sirnply replacing the trace 
by the relevant supertrace |2i l3i lj| i^ that the structure 
of the equations of motion and all other identities used 
are the same as for the np flavor case of refs. |ll[ Il2l | . 
However, the Cayley-Hamilton identity used there to re- 
duce the Lagrangians for the two and three flavour case 
is not valid in PQxPT. 

The last technical difficulty in PQxPT as compared to 
the mass calculations at two loops in xPT [lj,ll3 are the 
extra two-loop integrals needed since the propagators in 
PQxPT can have a double pole structure. These extra 
integrals can be evaluated using the same methods as was 
used for the integrals in ref. [lj|, provided that deriva- 
tives of those expressions w.r.t. the relevant masses in 
the propagators are taken. 



the integrals is the same as in ref. [lj| except that an 
extra integer argument now indicates the needed prop- 
agator structure. Index (1) corresponds to the case of 
single propagators only, as in ref. |l4l | , whereas index (2) 
indicates that the first propagator appears squared, in- 
dex (3) that the second propagator appears squared, and 
finally index (5) that the first and second propagators 
appear squared. Explicit expressions can be found in 
ref. [lj|, and by taking derivatives w.r.t. the masses of 
the expressions there. 
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FIG. 1: Feynman diagrams at 0{p^) or two-loop for the pseu- 
doscalar meson mass. Filled circles denote vertices of the C2 
Lagrangian, whereas open squares and shaded diamonds de- 
note vertices of the £4, and £e Lagrangians, respectively. 



THE PSEUDOSCALAR MESON MASS TO 0{p^) 

The mass of the pseudoscalar meson is obtained by 
means of dimensional regularization from the diagrams 
in Fig.[T]and those at ©(p^) and Oip'^). 

The result is expressed in terms of the valence quark 
mass, mqv, and the sea quark mass, mqs, via xi = 
2Bomqv, X4 = 2Bomqs and Xi4 = (xi + X4)/2. These 
correspond to the lowest order charged meson masses. 
Other parameters include the decay constant in the chi- 
ral limit (i^o): the quark condensate in the chiral limit, 
via {qq) = -BoFg, and the LECs of 0{p'^) and 0{p^), 
the LI and X[. 

The finite parts of the loop integrals in the expressions 
below are 



Mx) 
s(x,x;0) 

c{x,x,x;0) 



~7ri6xlog(x/M^), 
-^16 (1 + log(x/M')) , 
-^i6/(2x) , 



(1) 



where the subtraction scale dependence has been moved 
into the loop integrals. We also define Trie = l/(167r^). 
The finite two- loop sunset integrals H^ , H[ , 7?|^ that 
appear in the top right diagram of Fig. ^ may be eval- 
uated using the methods of [lj|. The notation used for 



We write the pseudoscalar meson mass in the form 



Mh = Xi 1 + ^^ 



c(6) 



loops 



p2 
-^0 



p4 



0(p8) 



(2) 



where we have separated the 0{p'^) and 0{p^) contri- 
butions. The two- loop part S'-^'> has been further split 
into the contributions from the chiral loops and from the 
0{p^) counterterms. The 0{p'^) mass shift is 

JW = +2Axi{2Ll-Ll) + 8xii2L;~Ll) 

-l/3i(xi)-l/3B(xi,Xi;0)[xi-X4], (3) 

which is in agreement with ref. JJ. The 0{p^) mass shift 
due to the Cq Lagrangian is 



r(6) _ 
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i2xlK{,- 

-48x1X4^20 

■16X?^2^3 + 
-AmxlKl,^ 



96xiX4i^[8-16x?^[9 
-48x^i^2^-144x4^22 
48x?if2"5+48[2xiX4 + xl]^; 
-32x?i^j9 + 96xiX4if: 



Mo- 



26 

(4) 



If one further defines 7ri2 — 7r^/12 -|- 1/2, then the mass 
shift from the remaining 0{p^) terms is 



loops 



TTJe [15/32 X1X4 - 3/32 xl + 73/64 xl] - ^\^ ^l^ [41/9 xiX^ + 10/9 x? + 17/4 x^] 
+ Trie [26/3 xiX4 - X? + 3x1] ^S + 4 ttis x? ^I + ^16 [2x? + 16x1] ^2 + ^i6 [17/3 xiX4 - 5/2 x? + 3/2 xl] ^5 
+ 3 ^16 A(xi4) X4 + 384 X1X4 ^4^5 - 1152 xl ^4^6 - 384 xiX4 L\Ll + 576 xl L? 

- 384 xiX4 LlLl - 128 x? LIL^ + 64 x? Lf - 8 ^(xi) [2xi - X4] L^ + 8 ^(xi) Xi L^ + 20 A(xi) Xi ^2 

- 8 i(xi) [2X1 - X4] LI + 16 A{xi) Xa L\ + 16 A{xi) [xi - 1/3x4] LI + 16 A(xi) [2xi - X4] LI 
+ 32 ^(xi) [xi - X4] l:, - 64/3 A{xi) Xi ^S - A{xif [13/72 - 5/18 xr'x4] 

+ ^(xi) 5(xi, Xi; 0) [77/36 Xl - 11/12 X4] - 2/9 A(xi) C(xi, Xi, Xi; 0) [xiX4 - xl] + 24 A(xi4) Xi4 ^S 
+ 60 A(xi4) X14 LI - 48 i(xi4) X14 LI + 96 ^(xm) Xi4 L^ - A(xi4)' [27/4 - xixr4'] 

- 2 A(xi4) B(xi, Xi; 0) X4 + 128 i(x4) X4 L^ + 32 A(x4) X4 L^ - 128 A(x4) X4 L\ 

? 

'21 L5 + S(xi,Xi;0) [32x1X4 -24x1] 

- B(xi, Xi; 0) [48x1X4 - 32x1] L^ - B(xi, Xi; 0) [32x1X4 - 16x? - 16x1] ^^ 

+ B(xi, Xi; 0) [64/3 xiX4 - 32x? + 16/3 xl] ig - S(xi, Xi; 0)' [1/3 xiX4 - 5/18 xl - 1/18 xl] 

+ B(xi, Xi;0)C'(xi,Xi,Xi;0) [2/9x1x1- 4/9 x?X4 + 2/9 X?]- 16 C(xi,Xi,Xi;0)[xixl-X?X4]i4 

- 16/3 (7(xi , Xl , Xl ; 0) [x?X4 - X?] LI + 32 C(xi , Xi , Xi ; 0) [xixl - X?X4] LI 

+ 32/3C(xi,Xi,Xi;0) [xiX4 -Xl] LI + 5/9 77^(1, xi,Xi, Xi; Xi) Xi -^^(l,Xi,Xi4,Xi4;Xi) [xi - 1/4X4] 

+ 2 ff^(l, Xi4, Xi4, X4; Xl) X4 - 4/9 H^{2, Xi, Xi, Xi; Xi) [XiX4 - Xl] 

-3/4iJ^(2,xi,Xi4,Xi4;Xi) [XiX4 -Xl] +2/9iI^(5,xi,Xi,Xi;Xi) [Xixl -2x?X4 + Xi] 

+ 4 -fff (3, xi4, Xl, X14; Xl) [xiX4 - Xl] + 3/4 i/fi(l, Xi, Xi4, X14; Xi) Xi 

+ 6ff|i(l,X4,Xi4,Xi4;Xi)xi +3/4i/|[(2,xi,Xi4,Xi4;Xi) [xiX4 -Xl]- (5) 

I 



128 ^(X4) X4 LI - A{xiY + 8/9 A(x4) 5(xi, Xi; 0) X4 + ^(xi, Xi; 0) [8x1X4 - 8x?] LI 
- 5(xi, Xi; 0) [8X1X4 -8x?]i3 + 5(Xi,Xi;0) [32x1X4- 24x1] L^-5(xi,Xi;0) [16X1X4 



56/3x?]iS 



As expected, the divergences in these expressions have 
cancelled. Note that in both eq. lO and eq. lO, the 
extra pole structure from the partially quenched nature 
disappears if we set xi = X4- The remaining B{xi, Xi! 0) 
terms at 0{p^) in that case are due to the fact that we 
have expressed the 0{p'^) result in terms of the lowest 
order masses rather than the full physical masses, which 
was called unrenormalized in ref. |l4l] . Consequently, 
they do not produce unphysical logarithms. 



DISCUSSION AND CONCLUSIONS 

In the long run, all input parameters should of course 
be determined by a fit of the formulas of PQxPT to 
lattice QCD data. At the present time, we can only 
present results using input parameters from the contin- 
uum work in xPT at the same order. We use as input 
parameters the fit to data presented in ref. [la], called 
fit 10, which had Fa = 87.7 MeV and fi = 770 MeV. Of 
the parameters which were not determined there, we set 



KI 



0. The last one cannot be deter- 



mined from experimental data and some recent results 
on L4 and Lg can be found in ref. |l6| . 

In Fig. 121 we show the hues in the Xi-X4 plane used 
in the following figures. The remaining figures show the 
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FIG. 2: The lines in the xi"X4 plane along which results in 
Figs. I3lf}l are shown. 



relative correction to the meson mass. In Fig. Owe plot 
M|g/xi - 1 to 0{p^), and in Fig. H to ©(/), in both 
cases with the LECs set to the fit 10 values of ref. |l3. 
Fig. also shows the result at 0{p^), but with the LECs 
set to zero. Note that the effect of the unphysical loga- 
rithms is clearly visible along the line labelled A which 
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has a constant sea quark mass. That the 0(j)^) con- 
tributions to the pseudoscalar meson masses are sizable 
was expected. Similar results for the case of xPT were 
obtained in refs. |l4l. Il,5l|. but cancellations with the con- 
tributions from the O(p^) Lagrangian might change this. 
In conclusion, we have calculated the off-diagonal pseu- 
doscalar meson mass at NNLO order in PQxPT, and 
have presented first results for the case of equal valence 
quark masses and equal sea quark masses. The 0(j}^) 
contributions arc sizable, and the effects of the loop con- 
tributions are definitely nonnegligible at presently used 
quark masses in lattice QCD calculations. 
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FIG. 5: The sum of the C'(p*) and £>(/) resuHs at ^ 
770 MeV with all LECs set to zero. 
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